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We investigate tiie two-dimensional classical dynamics of 
the scattering of point particles by two periodically oscillating 
disks. The dynamics exhibits regular and chaotic scattering 
properties, as a function of the initial conditions and parame- 
ter values of the system. The energy is not conserved since the 
particles can gain and loose energy from the collisions with 
the disks. We find that for incident particles whose velocity is 
on the order of the oscillating disk velocity, the energy of the 
exiting particles displays non-monotonic gaps of allowed ener- 
gies, and the distribution of exiting particle velocities shows 
significant fluctuations in the low energy regime. We also 
considered the case when the initial velocity distribution is 
Gaussian, and found that for high energies the exit velocity 
distribution is Gaussian with the same mean and variance. 
When the initial particle velocities are in the irregular regime 
the exit velocity distribution is Gaussian but with a smaller 
mean and variance. The latter result can be understood as 
an example of stochastic cooling. In the intermediate regime 
the exit velocity distribution differs significantly from Gaus- 
sian. A comparison of the results presented in this paper to 
previous chaotic static scattering problems is also discussed. 



I. INTRODUCTION 

The study of nonlinear systems capable of exhibiting 
chaotic behavior has been an intensive area of research in 
the last fifteen years. This research was initiated about 
a century ago by the work of Henri Poincare, who stud- 
ied the motion of three gravitationally interacting bod- 
ies. Most of the work done in this subject has focused on 
bounded systems. On the other hand, many experimen- 
tal techniques involve scattering processes. In contrast to 
bounded systems, where the particle's trajectories remain 
forever inside the range of interaction, in a scattering pro- 
cess an incoming particle feels the interaction potential 
only for a finite amaujit of time and eventually exits the 
interaction regionErO. In the general description of a 
scattering process, we have an input trajectory into a re- 
gion of nontrivial dynamics called the scattering region, 
and an output trajectory away from this region. We can 
think of the scattering process as a map that transforms 
an incoming trajectory into an outgoing one. Only rela- 
tively recently it has been realized that a scattering pro- 



cesses from a general scattering potential, often without 
a simple geometric symmetry, can have rather compli- 
cated dependencies between the incoming and outgoing 
trajectories. This means that, by very slightly changing 
the initial conditions that define the incoming trajectory, 
the outgoing one will have rather large fluctuations. The 
idea that chaotic scattering can play an important role 
in various problems in physics became wi|d|aly accepted 
after the seminal work of Petit and HenonQEl. 

Most previous chaotic scattering studies have assumed 
a stationary scattering region, i.e. fixed in time (for an 
exception sedl3). In this paper we present results from 
a dynamical study of the scattering of particles from a 
time-dependent oscillatory interaction potential, which 
consists of two circular disks that oscillate periodically in 
time. The static two-diskjiroblem was recently shown to 
be analytically integrableEiJ, (hereafter we call this work 
I). In this paper we build our nonequilibrium dynamical 
study based upon the results obtained in I. 

Our model can conceivably be produced in very low 
temperature experiments where a couple of circular quan- 
tum dots are generated by a gate voltage that can vary 
their radius periodically in time. Ballistic transport ex- 
periments in mesoscopic systems have raised the pos- 
sibility of directly studying chaotic billiards, where the 
addition of external fields can yield results that are ex- 
pected to account fox. certain aspects of unusual related 
experimental resultsEjo. Some of the transport results 
seen in experimeats are surmised to have classical related. 
explanationsEirEj. . The geometry of microjunctionslij 
and antidot-latticestSEj can be described by models that 
consist of circular scattering disks. For the above reasons, 
we will focus in this paper on the scattering of a parti- 
cle from two oscillating classical hard-disk billiards. Here 
we concentrate on the classical dynamics of this model 
and leave the very interesting quantum case for a future 
study. 

The outline of the paper is the following: In Section II 
we introduce the model considered in this paper, together 
with its main physical properties. In section III we derive 
a scattering map associated with our problem. In section 
IV we present and discuss the bulk of our results. Finally, 
in section V we provide a short summary of the results 
and the perspectives for the future. 
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II. DEFINITION OF THE PROBLEM. 

We consider the motion of a unit mass particle re- 
stricted to move on the plane. The particle elastically 
collides with two hard disks that oscillate periodically in 
time. The initial velocity of the particle changes as a 
function of time due to the energy exchange after each 
collision with the disks. As we discuss below, depending 
on the initial conditions the particle will spend a cer- 
tain amount of dwell time in the interaction region, after 
which time it will exit upwards or downwards. It is the 
complexity of this motion that we will carefully describe 
below. _ 

Here we will follow the approach presented in lO, in- 
cluding its notation. The reader should check this refer- 
ence for further details on the formulation of the static 
problem. In Fig. we show the two disks on the 
plane. The radius of both disks are normalized to one. 
Their centers are separated by a time-dependent distance 
R{t) > 2. One convenient way to study this problem, as 
pointed out in I, is by replacing the system by one disk 
and one rigid wall placed at the symmetry axis of the two- 
disk problem. This is the representation of the model we 
study in this paper. 



A. Two-disk oscillating model. 

The model we consider here is in some sense the 
scattering two-dimensional extenskin of the well studied 
bound Fermi acceleration modelEHo. This model is de- 
fined by a free particle inside a rigid one-dimensional box, 
with one wall fixed and the other one periodically oscil- 
lating in time. The Fermi model was one of the first two 
degrees of freedom problems studied, which exhibited a 
transition from regular to chaotic behavior as a function 
of the oscillating wall motion. For a linear saw-tooth 
time-dependent wall oscillation, the particle dynamics is 
regular. Having a linear time-dependence implies a con- 
stant oscillating wall velocity. When the oscillation is 
nonlinear in time, there is acceleration in the wall mo- 
tion and one can then have non-trivial dynamics, with a 
transition between regular to fully chaotic behavior. In 
this paper, without loss of generality, we consider the 
simplest nonlinear piece-wise quadratic time-dependent 
disk oscillation shown in Fig. (|2|). In this case we repre- 
sent the motion of the disk center by 



integer m each continuous piece of the disk oscillation. 
The time t and m are then related by the expression 



X{t) ^ At^ + Bt + C. 



(1) 



Here the constants A, B and C are fixed for a half period, 
and they have different values for different half periods. 
The motion of the disk center (the other disk is the mir- 
ror image of this one) is given by Eq. modulus the 
oscillation disk period T. 

In our analysis, for calculational convenience, we chose 
to treat the problem in the following way. We label by the 



2{t + <t>o/uj) 
T 



(2) 



where | ] denotes the nearest lower integer. The m pa- 
rameter will have a fixed value for time t in the time 
interval if)T<t< (^)T. Here uj = 2tt/T is the os- 
cillation frequency, and $o is the initial oscillation phase. 
The specific expressions that define the parameters A, B 
and C are given in Appendix A. 

We have now defined the time dependence of the oscil- 
lating disk. Next we use the relevant results given in I, 
noting that the incidence-reflexion symmetry in our case 
is changed by the oscillation of the disk. 



1. Collision Time 

We start by calculating the time elapsed between two 
successive collisions of the particle with the disk. We 
need this time to calculate the new velocity vector, by 
means of a velocity transformation to the system where 
the disk is at rest. We deduce from Fig. that the 
position of the colliding particle is given by 



p{tn+l) = Pn+l ^ Pnv + ^ 



(3) 



where tn is the previous collision time, and the subindex 
V denotes a specular variable. To clarify the meaning of 
specular, consider for example, the one associated with 
Pn{t) = (a;„,y„) which gives 

Pnv — {p^nviDnv) — ( •^n-iUn^ — ( -^n H" COS ^tj,, sin ^j^,), 

(4) 



and the velocity 

= {vn cos(6'„ - 0„), w„ sin(6l„ - 0„)). (5) 

The new collision point Pn+i = (a^n-i-i, J/n-t-i), at the new 
collision time must lie on the circumference given 

by the equation 



(a;„+i - X{tn+\)f + y\+x = 1- 



(6) 



Evaluating Eq. (|l|) at i = tn+i, and substituting it in 
Eq. (^), we get the quartic equation for tn+i 

aitn+i + aai^+i + a2<fi+i + ai^n+i + ao = 0, (7) 

where the expression for the parameters aq, ai, 02, 03 
and 04, are explicitly given in Eq. (A. 2). 

We can get tn+i as a function of t„ using Eq. and 
Eq. (A. 2). Once we know the collision point on the 
disk can be determined from Eq. (P) as 



2 



(8) 

and then, using Eq. (O), the disk will be located at 



Using Eq. (JlSj) in the last equation we get 



1+1 



(9) 



2. Disk velocity Map. 

To calculate the velocity of the disk, Vn+i, at the new 
collision time, we take the time-derivative of Eq. (|l|) that 
gives 

V{t) = {X{t),0) = {2At + B,0), (10) 
and consequently 

Vn+l = V{tn+l) = (2it„+i +B,0),= (K+1,0) (11) 

which is fully determined since is known from Eq. 
(^. To determine the velocity of the particle, Vn+i, we 
introduce the relative particle velocity (see Fig. (||)) with 
respect to the disk as 



Then 



where 



and 



u = V — V . 



Unv A fn+1 = Pfc, 



P = {{Xn+1 - Xn+l)Uny - Vn+lUnx), 



(12) 
(13) 
(14) 



with the normal unit vector k \\ z, as seen in the figure. 
One can also show that 

Un+l A f„+i = {Un+ixUn+l " Un+ly{Xn+l - X„+i))fc. 

(16) 

From these two equations we get 

Un+ly = P - QUn+lx, (17) 

with 



P = 



Xn+l — Xn+1 



and Q = 



Vn+l 



(18) 



We also have the conservation of the velocity magnitude, 
which in the coordinate frame where the disk is at rest is 
given by 



2 I 2 



n+lx ~^ "n+ly ~ 



(19) 



Un+lx 



PQ-^P^Q^-il + Q^){P^ 

1 + Q2 



(20) 



with an appropriately chosen " — " sign in the square 
root. We now go back to Eq. ( p^ to find Un+iy This 
allow us to obtain the velocity of the particle after the 
collision through the expressions 

^n+lx — '^n+lx + Vn+1 and Un+ly = Un+ly (21) 



III. THE SCATTERING MAP. 

Following the notation of I, we can get the scattering 
map associated with this dynamical system. Since the 
derivation of our map is completely analogous to the one 
given there, we can directly write down the final expres- 
sions stating a few differences proper to our problem. 
The map derived in I is 



^1 = sm 



^n+l 



-(sin(/)„ + i?sin(6'„ - </)„)) 



(22) 



0„+i = sin ^ {sin 9n + Xsin{9n - ifin))- (23) 
In our case R ~ Xn + ^n+i and 
A = Rcos{9n 



^n V^n I 



[cos(/)„ - i?cos(6in - 0n)]^ - + 2i?cos6i„. (24) 



The initial conditions {9o,(f)o) in this map are obtained 
from the parameters set at time t — 0. If we take the 
initial particle position as (xq^i/q), the angle a between 
the initial velocity Vq and the horizontal gives 



sin ^ 



{{Xq — Xq) sin a + j/g cos a) 



(25) 



6*0 = sin {uq + Ag sin a). (26) 
In Eq. (^5|), Vq is given by vq = uq + Vq, and uq by 
RS - Jr^S^ - (1 + 5-2) (i?2 _ uD 



UOx 



1 + 



UOy ~ R~ Suox, 



(27) 



(28) 



where 



R ^ and S ^ , (29) 



Xq - Xo 



Xq - Xq' 
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R= (Xq- XQ)uQy + yoUQ^ with u^^ = Vq^ - Vq 

(30) 

The results derived in this section used a polar coordi- 
nates representation that has several advantages for the 
geometric analysis described here. When iterating the 
map numerically, however, the polar coordinate repre- 
sentation is somewhat cumbersome and for that reason 
we found it more convenient to carry out the iterations 
in Cartesian coordinates. This is what we did to obtain 
the results described in the next Section. 



IV. RESULTS 

In this Section we discuss the bulk of our numerical 
results. We provide typical results for a regime of in- 
teresting physical parameters. To check our analysis, we 
looked at the Fermi acceleration limit of our problem, 
which corresponds to having the two disks quite close to 
each other and with the particle initial conditions along 
the disks axis, so that the particle does not notice the 
disks curvature. We reproduced the Fermi accelerator 
model results by choosing the parameters for the equi- 
librium position of the disk center, X^, the amplitude of 
oscillation. A, the time oscillation period T, and the free 
space distance betwccjp-thc wall and the disk, close to the 
values given in Refs.EZllHj. The phase space plots obtained 
correspond well to the known Fermi accelerator results, 
i.e. chaotic behavior for low velocities, and several sets 
of resonant islands for higher velocities. 

After this test we proceed to chose the separation be- 
tween the wall and disk large enough so that the particle 
dynamics sensed the curvature of the disks. Of course, 
if the separation distance is too large the particle will 
hardly collide with the disk and the dynamics becomes 
trivial. The interesting parameter ranges are the ones 
which allow a large number of particle collisions with the 
oscillating disk and the wall. A typical phase space plot is 
shown in Fig. (^, for several particle initial conditions. 
The parameters considered satisfy the necessary condi- 
tion to have a large number of particle collisions with the 
oscillating disk. In the units where the radius of the disk 
is one, we took the parameters: (Xe,0) = (1.000097,0), 
A = 1.6 X 10-^, T = 7.6 X 10"^ $o = 0, and the accel- 
eration parameter A = 4716.085. We considered a set of 
2000 particle initial conditions, that we may as well call 
a beam of 2000 particles, each one sent from the origin 
into the scattering region with an angle a = 6 x 10~^ ra- 
dians with respect to the x-axis. We varied the velocities 
of the particles between 0.1-5.0, chosen from a uniform 
random distribution. 

In Fig. (||) we show the delay, or dwell time, t^, as 
a function of the initial energy (velocity) of the incident 
particles. For low energies we observe a very irregular 
behavior in Td- In fact, this behavior is rather close to a 
fractal, as can be seen from zooming in a given interval 



•8f^ener^8s,^^Sjiown in the left inset. We also note that 
for initial velocities larger that 5.5, there is a mixture of 
regular and irregular zones. When we amplify one of the 
irregular regions, we can again see the fractal character 
of the results (see inset on the right). For larger ener- 
gies than the ones shown here, we found that tends 
to a constant value. This is what we would expect for 
large energies since the particle essentially sees the disk 
as stationary. 

It is interesting to use the data of Fig. (||) to con- 
struct the histogram of dwell times shown in Fig. (|^). 
The main histogram has two representative contribu- 
tions. One comes from the irregular zone and the other 
one from the semiregular component, as it is shown in 
the two insets in the figure. The upper inset corresponds 
to the irregular region and the lower one to the semireg- 
ular zone. The main histogram shows one peak close to 
a Td of about 100 and the other one close to 200, which 
correspond to the peaks seen in the insets. The bin size 
used in all the histograms shown are around 3% of the 
full range. In Fig. (|^) we show the number of collisions 
with the disk versus the incoming velocity. The general 
behavior is very similar to that in Fig. (|5|), however, the 
explicit relationship is complicated. This can also be seen 
from Fig. (||), which we shall discuss in the next para- 
graph. In Fig. (|^) we note, in particular, the irregular 
behavior in the same regions of incident velocities. As we 
increase the initial energy, the number of collisions reach 
a plateau, which is because the disk appears to be at rest. 
In the inset we used the same number of particles as used 
in the main figure, but the range of velocities is smaller, 
so as to allow us to see in more detail the results. 

In Fig. (^ we display the dwell time vs the number 
of collisions, N, to see if there is a simple relation be- 
tween them. The input velocity appears here only as a 
hidden variable. For example, when the range of veloci- 
ties in the inset is in the irregular region, there is a wide 
spread in the location of the resulting points. If we also 
allow initial velocities from the semiregular region, as it 
happens in the main figure, then the data points are lo- 
calized around a specific zone that is darker in the figure. 
This is consistent with an essential independence of these 
variables when the initial velocities are outside of the ir- 
regular region. The number of collision data points in 
the main figure and in the inset are equal. 

Next we discuss the relevant scattering variables of the 
problem. In Fig. (^ we show the irregular behavior of 
the exit angle as a function of the normalized initial par- 
ticle velocity. The low energy particles, with velocities 
less than 5.3, get an irregular exit angle in a wider range 
of values. When we plot the distribution of these exit 
angles, we find a wide pattern centered around an an- 
gle of about 0.45 radians. This can be seen in the inset 
of this Fig. (H). Particles with incident velocities larger 
than 5.3, that show semiregular behavior, also contribute 
strongly to this peak. The corresponding histograms for 
these regions are displayed in Fig. (p^), with the left his- 
togram associated with the irregular region, and the right 
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one with the semiregular region. Both histograms show 
a peak for an exit angle around 0.45 radians. For low 
energy incident particles we get a wider range of output 
angles. 

In Fig. ([l^) wc show the exit velocity as function of 
input velocity. The region with input velocity less than 
5.3 is quite irregular, and it is consistent with the previ- 
ous figures. When we increase the input velocity, the exit 
velocity grows and the fluctuations are about a line with 
slope of almost 45°. In this figure it is noticeable that 
there are big jumps for input velocities between 16.7 to 
19. The left inset is an amplification for low input veloc- 
ities. The right inset shows the corresponding histogram 
for the exit velocities. Here we notice that there are iso- 
lated peaks or gaps in this distribution. 

The histograms shown in Fig. (|l^) are directly related 
to the data shown in Fig. (|ll| ) . The left histogram corre- 
sponds to the data given in the inset of Fig. ([ll|), while 
the right one corresponds to the main plot. In both ana- 
lyzed histograms the fluctuations are about the line with 
slope 7r/4, with data points on this line labeled by the 
variable v. The variables V^ax and Vout are the max- 
imum disk velocity and the real output velocity of the 
particles, respectively. One important feature of these 
histograms is that they appear to be directly related to 
the isolated peaks mentioned before. Here the effect is 
more prominent for the irregular region of input veloc- 
ities. The effect remains, even if we make the velocity 
equal, but with different initial phases. This is equiva- 
lent to carry out a phase average in the interval (0, 27r). 
The gaps in the output velocities are also gaps in out- 
put energy. These results indicate that there are output 
energy regions that the particle can not explore, leading 
to forbidden energy regions. In the inset at the top of 
Fig. (pl|), we note that the exit velocities have a peak 
when the input velocities are close to 19. For the range 
of velocities between 23 — 60 the energy gaps were not 
seen. This is why the gaps are wider at the bottom in the 
inset of Fig. (|ll|). If we increase the range of input ve- 
locities, as in the main figure, then there appear narrow 
gaps related to different velocity contributions. 

We have also carried out a basic fractal analysis of 
a 10,000 particle system. The idea was to extend the 
analysis of Ref.cl to two-dimensions. We determined the 
plane boundary of initial conditions {xq^Uq), which sep- 
arates the particles into the ones that go upwards from 
the ones that go downwards. We plotted a figure with 
black squares representing the initial conditions of par- 
ticles which exit upwards, and empty squares belonging 
to the ones that go downwards. We obtained 1.86 as the 
fractal dimension. We do not show these results since 
they are typical of chaotic scattering problems. We car- 
ried out this quantitative analysis to make sure that all 
the qualitative generic properties of a chaotic scattering 
system applied. Although all the results are quantita- 
tively different, as one should expect, we did not find 
a significant change in the general qualitative behavior 
described above. 



Finally, we note that the model we are considering here 
does not conserve energy, and we are also interested in 
understanding how energy is added or subtracted from 
the disk to the colliding particles. One possibility is to 
take the initial velocities distributed by a Gaussian func- 
tion, just as in the classical statistical mechanics Maxwell 
velocity distribution. We chose then a beam of particles 
with this velocity distribution with a given standard de- 
viation (T, or inverse temperature. Then we studied the 
evolution of the distribution of exit velocities. We did 
the analysis at low, intermediate and high beam ener- 
gies. The results are shown in Fig. (|l^). In all of these 
figures, the continuous line curve represents the Gaussian 
distribution fit to the beam of incident particle velocities. 
The histogram is the distribution of exit velocities after 
scattering. We note that the Gaussian distribution is 
maintained only for high energies (right-bottom figure), 
but for low or medium energies the exit velocities can be 
not be fitted to a simple Gaussian. At low energies, when 
the incident velocities are in the irregular region (left- 
top figure), however, most particles concentrate about 
a Gaussian-like distribution, with smaller mean and a. 
These latter results indicate that the beam losses energy 
and that it has some kind of stochastic cooling. 

V. CONCLUSIONS. 

In the present paper we considered the complex dy- 
namics of a particle that scatters from two periodically 
oscillating disks with a variety of initial conditions. We 
found that the dynamics has regular and irregular be- 
havior that we analyzed in some detail. This model is in 
a sease a dynamical extension of the well known Lorentz 
gasEj. Although the model studied here is perhaps the 
first dynamic chaotic scattering analysis, several of the 
results we described are similar to those found in chaotic 
static scatters. There are, however, some important un- 
expected differences in the results obtained. Among the 
most interesting and surprising results presented in this 
paper are the energy gaps found in the exit energy. This 
result indicates that there is an important energy ab- 
sorbing mechanisms, directly related to the nature of the 
classical dynamics of the problem. We found that the 
markedly irregular dynamics appears when the particles 
have velocities on the order of the disk velocities. It is 
within this energy range that the energy gaps appear. 
For larger particle energy the dynamics simplifies, for 
the oscillating disks appear as if they were at rest. 

Another important difference from the dynamics of 
chaotic static scatterers, has to do with the energy gained 
or lost by the beam of particles. As a test, we took a 
Maxwell-like distribution of initial velocities and found 
that only in the large energy regime, where the disk is 
essentially seen at rest by the particles, the Gaussian dis- 
tribution of exit velocities is preserved. Otherwise, there 
are important changes in the exit velocity distributions 
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for low velocities. 

In this paper we considered that the disks do not ab- 
sorb energy from the colliding particles. Including energy 
gained or lost by the disks is necessary to mimic the ef- 
fect of temperature in the model. We intend to include 
this effect in the model elsewhere. The very interesting 
questions raised by the quantum mechanical treatment 
of the model introduced here are left for the future. 
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APPENDIX: A 

In this appendix we write down the explicit expres- 
sions for the parameters defined in the main body of the 
text. In order to obtain the values of the parameters A, 
B and C of section II, we use Fig. (||), and ask that the 
parabohc curve crosses trough the points (^T, + A) 
and (^^^T, Xe — A), where Xe is the equilibrium po- 
sition of the center of the disk and A is the oscillation 
amplitude. When this is done, we obtain the following 
expressions: 



The work by AA was partially supported by DGAPA- 
UNAM grant IN-105595, and the one by JVJ was sup- 
ported in part by NSF grant DMR-95-21845. 



A ~ free parameter , 

OA ( 4\ 

B = i-ir— + - (2<I>o - ^(1 + 2m)) 

TT \ U 



C = Xe + (-1)"+M(1 + 2m) + (-1)' 



,2A<S>o 



A 



— ($0 + W7r2(l + m) - 7r$o(l + 2m)) 



(Al) 

A is the curvature in the saw tooth, which is associated 
to the disk acceleration. 

The derivation of the parameters ao, ai, a2, and 
04 appearing in Eq. (|^) is a straightforward and here we 
just cite the results. We evaluate Eq. (|l|) at t = tn+i and 
then substitute the result into Eq. (||). Then we obtain 
Eq. (Q) with the following coefficients. 

ao = a;„ + Vn^n + C* — 2XnVnxtn + 2x„C' 

ai = -2vljn + '2BC' + 2XnVnx + 2XnB 

'^^nx^^n ~t~ 2f^^3;C* 2Vy^ytji -|- 2yjiVjiy^ 

02 = + B^ + 2AC + 2x^A - 2vnxAt^ + 2f„xi?, 
as = 2AB + 2vnxA, 
04 = A^. 



{A.2) 
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FIG. 1. This figure defines the two-disk periodically oscil- 
lating model studied in this paper. The model is replaced and 
studied by the one oscillating disk and a fixed wall. The fig- 
ure shows a dotted line disk with center at X„ = and 
a continuous line circle at the right of the wall with center at 
Xn+i = X{tn+i)- The disk on the left of the wall represents 
the image of the dotted line disk. Variables with a subindex 
n are evaluated at time tn- 



FIG. 2. Here X{t) denotes the oscillating wall 

model about the equilibrium position Xe in the range 
[Xe — A, Xe + A] , with amplitude A. The parameter rn labels 
the oscillation segment with fixed value in the time interval 
i^T, -^^^T). The relation between time t and m is given in 
Eq. (2). The figure is drawn for $o = 0. 



FIG. 3. This figure is used in the derivation of the velocity 
map. The relative velocity before and after the collision has 
the same angle with respect to the normal to the disk, but 
the velocity itself has two different angles. See text for the 
definition of the variables in this figure. 



FIG. 7. Here we show the dependence between N and 
Vin/Vmax- The main plot and the inset show the same general 
behavior as that of Fig. (5), with the irregular behavior in 
the same range of input velocities. For largo velocities, A'^ is 
basically constant. This means that the disk motion has little 
effect in the particle velocities. Each picture was done with 
2000 particles. 



FIG. 8. Here we show a plot of ra vs N. We note that 
there is no clear relationship between these two variables. The 
general behavior of both plots arc similar to the ones shown 
in Figs. (5) and (7). Note that for low Vin (irregular region), 
their dependence is irregular. The number of points taken 
was 2000. We see that when the input velocity is outside the 
irregular region, all the particles concentrate in a black zone 
in the figure. This moans the almost independence of ra and 
N with respect to the input velocity. 



FIG. 4. Phase space results for a limit close to the 
Fermi acceleration one-dimensional model. Here Vx is the 
x component of the particle velocity, normalized by the disk 
velocity Vmax- The particle will eventually sense the disk 
two-dimensional curvature, and the resonances structure, as 
compared to the Fermi model, will change. The plot has 
74238 points obtained from 17 particles with different initial 
velocities. 



FIG. 5. Dwell time Td as a function of the normalized 
incident velocity Vin. In inset a we show amplified results 
about Vin = 2.375. We see that the low energy particles have 
an irregular behavior with a fractal-like character. Velocities 
higher than 5.5 have a mixture of regular and irregular be- 
havior. In inset 6 wo show us a further amplification around 
velocities close to 7.581. Each picture is drawn from results 
of 2000 initial conditions. 



FIG. 6. Dwell time histogram N{Td/T) for data like the 
one shown in Fig. (5), except that here we used 6000 particles 
instead of the 2000 used in Fig. (5). The two noticeable 
peaks are due to the contributions from the irregular and 
semiregular zones, which correspond to the data shown in 
insets a and 6, respectively. The number of particles used to 
obtain the histograms was 2000. The bin size in the three 
histograms is around 3% of the full range in each plot. 



FIG. 9. Exit angle Qout as a function of Vm/Vmax- The 
low velocity region again displays an irregular behavior for 
the exit angle. The inset in the figure for the preferred exit 
angle has a peak around 0.45 radians. 



FIG. 10. Histograms of the data N{aout) shown in the 
inset of Fig. (9) that separate irregular from semiregular re- 
gions. Both histograms have a peak close to Qout = 0.45 
radians. The irregular region histogram is shown in a. It 
covers a wider range of output velocities than the semiregu- 
lar region shown in histogram b. The bin sizes taken for the 
histograms are 3% of the full range in Oout- 



FIG. 11. In the main figure we show the exit velocity vs 
input velocity for 6000 particles for a wide range of initial ve- 
locities. Inset a has the same coordinates, but for 2000 input 
velocities. In inset a the irregular region is more detailed, and 
between 0-15 the coarse averaged slope is small, while outside 
this range the averaged slope is close to 7r/4. We use the 
latter result in the histogram of Fig. (12). The exit velocity 
histogram is shown in inset b. Notice the isolated peaks in 
the distribution. 
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FIG. 12. Histogram of the N{vout — v)/Vmax data, a 
corresponds to the data given in inset of Fig. (11) (a), and 
b to the data shown in the main plot. The fluctuations arc 
analyzed around the line with slope n/4, and are labeled by 
the variable v. Here Vout is the real output velocity of the 
particles. We surmise that some of the isolated peaks present 
may be connected with the exit energy gap regions. 



FIG. 13. Each of the four figures was produced from a 

beam of particles with a Gaussian initial velocity distribution, 
denoted by a continuous line in the flgures. The histograms 
of N{vi„) or N{vout were produced from the exit velocities 
after scattering the oscillating disks. Each figure has a dif- 
ferent mean and standard deviation. In figure d we see that 
the Gaussian distribution is preserved but only for high en- 
ergies. In the irregular exit velocity region (a), the exit dis- 
tribution is concentrated about a Gaussian-like distribution 
with smaller mean and a. This result can be interpreted as 
some kind of stochastic cooling mean and a. This result can 
be interpreted as some kind of stochastic cooling Figures (b) 
and (c) correspond to intermediate energy regimes where the 
exit distribution is not Gaussian. 
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